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We use an optimal estimator to study the variance of the WMAP 9 CMB field at low res- 
olution, in both temperature and polarization. Employing realistic Monte Carlo simulation, 
^ we find statistically significant deviations from the ACDM model in several sky cuts for the 

Qh temperature field. For the largest mask considered in this analysis, which leaves 46% of the 

O sky, the CMB sky and ACDM are incompatible at > 99.94% C.L. at large angles (> 5°). We 

c/2 find instead no anomaly in polarization. As a byproduct of our analysis, we present new, 

a 

'~~' optimal estimates of the WMAP 9 CMB angular power spectra from the WMAP 9 year data 



at low resolution. 



INTRODUCTION 



The cosmic microwave background (CMB) anisotropy field can be used to probe cosmology 
►>■ to high precision, as shown by the Wilkinson Microwave Anisotropy Probe (WMAP) 9 years 

rS results |13j and by the very recent Planck cosmological results (see [21] and references therein). 

CMB data have given a significant contribution in setting up the A cold dark matter (ACDM) 
cosmological concordance model. The latter establishes a set of basic quantities for which CMB 
observations and other cosmological and astrophysical data-sets agreq^ spatial curvature close 
to zero; ~ 68.5% of the cosmic density in the form of Dark Energy; ~ 26.5% in cold dark matter; 
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^ See, however, |24| tor a tension concerning the Q.m estimate from Planck CMB and galaxy clusters data. 



~ 5% in baryonic matter; and non perfectly scale invariant adiabatic, primordial perturbations 
compatible with Gaussianity [23j. The largest angular scales of the CMB anisotropics have 
remained largely primordial, having evolved only in amplitude since the end of the inflation. 
They can be used to probe the physics of the very early Universe. Accurate estimates of their 
amplitude is a crucial to extract information about such primordial phases |25j . 

In this paper we provide new, temperature and polarization, estimates of the CMB Angular 
Power Spectra (APS) from the WMAP 9 year low resolution data using a Quadratic Maximum 
Likelihood (QML) method [31]. We have used BolPol, an efficient, fully polarized, implementa- 
tion of the QML method. BolPol has been already applied to WMAP 5 [lU] and WMAP 7 [II] 
low resolution data. At large angular scales, QML estimates are superior to the pseudo spec- 
trum techniques normally used in the CMB world |9l [T^ [271 EH] because of their smaller intrinsic 
variance |T7]. The WMAP 9 analysis has been derived using a suboptimal implementation of 
the QML method that extends to high multipole [Ij. BolPol is fully optimal, but computational 
considerations limit its application to smaller multipoles. This is fully adequate for the low res- 
olution analysis presented here. For other applications of our QML implementation see |191l20j. 
where it is used to constrain the CMB hemispherical power asymmetry in WMAP 5 and 7 data 
respectively, and [12] , where the low resolution spectrum of the cosmological birefringence angle 
has been derived for the first timq3 

In this paper, we provide new APS estimates from WMAP 9 year low resolution data. At 
the same time, we check for their stability against Galactic masking. Specifically, we look at the 
variance of the temperature and polarization maps. The variance of the CMB temperature field 
has been studied in previous works (see [3j for WMAP 5 and 7 year data and [18j for WMAP 
3 year data) which, however, focused on the small scale pattern. Here we focus on the lowest 
multipoles of WMAP 9 year data, we adopt a new estimator (based on BolPol) and consider 
also polarization. 

The paper is organized as follows. In Section IIIj we briefiy describe the QML method and 



our implementation, BolPol. Section III is devoted to the description of the WMAP 9 year low 



resolution data set under consideration. APS estimates of this data set are given in Section IV 
where the consistency of all the spectra with what expected in the ACDM model is evaluated 
through a x^ analysis. In the same Section the stability of the obtained estimates is studied 
against the Galactic masking. The analysis of the variance is presented in Section IVl Our 



conclusions are drawn in Section VI where some speculations about a possible origin for this 



2 



See also [^, where BolISW , a code largely based on BolPol, has been used to detect the ISW effect 



effect are set forth. 

II. ANGULAR POWER SPECTRA ESTIMATOR 

In order to evaluate the APS we adopt the QML estimator, introduced in [30J and extended 
to polarization in [31J. In this section we describe the essence of such a method. 

Given a map in temperature and polarization x = (T,Q,U), the QML provides estimates 
Cf - with X being one of TT, EE, TE, BB, TB, EB - of the APS as: 

Cf = E (P-')V [x*E^,x - tr(NE^O] , (1) 

£' ,X' 

where the Fyv' is the Fisher matrix, defined as 



F^'x' = \tr 



_i dC . dC 



C"^— -^C 



9C/^ dc'f' 



and the E C matrix is given by 



r.£ ^ i~>-l ^^ i~>-l 



with C = S{Cf) + N being the global covariance matrix (signal plus noise contribution). 

Although an initial assumption for a fiducial power spectrum Cf is needed, the QML method 
provides unbiased estimates of the power spectrum contained in the map regardless of the initial 
guess, 

(Cf) = Cf , (4) 

where the average is taken over the ensemble of realizations (or, in a practical test, over Monte 
Carlo realizations ex-tracted from Cf). On the other hand, the covariance matrix associated 
to the estimates, 

{ACfACf) = {F-')V , (5) 

does depend on the initial assumption for Cf: the closer the guess to the true power spectrum 
is, the closer are the error bars to minimum variance. According to the Cramer-Rao inequality, 
which sets a limit to the accuracy of an estimator, Eq. ([5]) tells us that the QML has the 
smallest error bars. The QML is then an 'optimal' estimator. 

This method has been implemented in a F90 code, named BolPol. Further details can be 
found in ITOl. 



III. DATA SET 

In this Section we describe the WMAP 9 year data set that we have considered. We use the 
temperature ILC map smoothed at 9.1285 degrees and reconstructed at HealPi5pl[H] resolution 
^side = 16, the foreground cleaned low resolution maps and the noise covariance matrix in 
(Q, U) publicly available at the LAMBDA website^ for the frequency channels Ka, Q and W as 
considered by [IJ and [16j for the low I analysis. These frequency channels have been co-added 
as follows dH US] 

mtot = Ctot{C]^lrnKa + Cq^^q + Cy^mv) , (6) 

where m,j, Cj are the polarization maps and covariances (for i = Ka, Q and V) and 

^tot = C'xa + ^Q +Cy ■ (7) 

This polarization data set has been extended to temperature considering the ILC map. We have 
added to the temperature map a random noise realization with variance of IfiK^ as suggested 
in [7J. Consistently, the noise covariance matrix for TT is taken to be diagonal with variance 
equal to IfiK^. 

IV. ESTIMATES 

Taking into account the kq85 mask for the temperature mapjj and the P06 mask for polar- 
ization maps we estimate all the 6 CMB spectra of WMAP 9 year data with BolPol, see red 
symbols of Fig. [T} In the same Figure we also show the averages of 10000 "CMB plus noise" 
Monte-Carlo (MC) simulations (see black symbols) where the CMB is randomly extracted from 
a ACDM model, drawn as dotted line, and the noise is obtained through a Cholesky decom- 
position of the noise covariance matrix. Consistently with WMAP data, also the temperature 
simulated maps have been smoothed at that scale, i.e. 9.1285 degrees. 

A. x^ consistency with ACDM model 

In order to evaluate the consistency with the ACDM model (identified in the following with 

a set of Cj for each spectrum) we adopt a frequentist approach, see Fig. ^ For each spectrum 

^ http://healpix.jpl.nasa.gov/ 

■* http://lambda.gsfc.nasa.gov/ 

^ Best-fit monopole and dipole have been subtracted from the observed ILC map through the HealPix routine 

REMOVE-DIPOLE [H]. 
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FIG. 1. WMAP 9 year CMB spectra at low resolution. ILC in Temperature, masked with kq85 mask 
and the co-addition of Ka, Q and V channel in Polarization masked with the P06 mask. From top to 
bottom, TT, EE, TE, BB, TB and EB spectrum. For each spectrum we show i{£ + l)Ce/2n (given in 
//K^) versus the multipole order £. Red estimates are for the WMAP 9 year maps. Black estimates are 
the averages of 10000 Monte-Carlo simulations where the CMB is randomly extracted from a ACDM 
model (shown as dotted lines) and the noise is obtained through a Cholesky decomposition of the noise 
covariance matrix. See also the text. 
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FIG. 2. x^ test for consistency with ACDM model. Each panel shows the counts versus x^- The 
histograms are built through Monte Carlo simulations. The red lines are the x^ values we obtain for the 
WMAP 9 year low resolution data. See also the text. 



(TT, EE, TE, BB, TB and EB) we plot the histogram of the xf values defined as 



Xi 



£,^'=2,32 



th\ 



(8) 



where i = 1, 10000 runs over the MC realizations and where the covariance matrix M^' is 
defined as 

Mw = ((Q,, - Cf){Ci,^i - Cf)) , (9) 

with (• • •) standing for the ensamble average. These histograms will be compared with the 
XwMAPg obtained with the WMAP 9 year data for each spectrum, defined as 

XwMApg = X! (C(?,W9 - Cl )M^g, (Q',w9 -Ce') , (10) 

^,£'=2,32 

where C^^wg are the WMAP 9 year spectra. Fig. [2] shows that the WMAP 9 year low resolution 
data appear to be well consistent with the ACDM model defined by the WMAP 9 fiducial 
model, represented by the dotted line shown in Fig. [T} See Table |l] for a quantitative analysis. 



TABLE I. Percentage to obtain a x^ value from a random extraction of ACDM model smaller than what 
observed by WMAP 9. 

Spectrum Percentage 



TT 


59.26 


EE 


88.99 


TE 


59.25 


BB 


46.13 


TB 


12.22 


EB 


85.73 



TABLE IL Sky fraction of the sky coverage of the considered masks. See also Fig. [3] 



Case Extension 


Sky fraction Sky fraction 


wrt kq85 and P06 


in T 


in P 


a +0 deg 


0.78 


0.73 


b +4 deg 


0.68 


0.65 


+8 deg 


0.56 


0.53 


d +12 deg 


0.46 


0.43 



B. Stability of the results versus Galactic masking 

In this Subsection we test the stability of the WMAP 9 year data estimates considering var- 
ious Galactic masks, see Fig. [3j More specifically, we extend the edges of the kq85 temperature 
mask and the P06 polarization mask by 4, 8 and 12 degrees. In Table |ll] we provide the sky 
fraction for each considered case and for reference we rename with the letters "a" , "b" , "c" and 
"d" the four cases we take into account. Case "a" is what is already considered in Subsection 

In Fig. [4] we show the six spectra for all the four considered cases. The overall consistency 
is good. However, the lowest multipoles of the TT spectrum show a clear trend: increasing the 
temperature mask, we obtain a lowering of the power at the largest scales from ^ = 2 to ^ ~ 7. 
Fig. [5] shows a zoom on the TT panel of Fig. |4] for the lowest multipoles, where the estimates 
are now joined with broken lines in order to make this behavior clearer. 




FIG. 3. Galactic masks. Upper panel: Masks we adopt in Temperature. Lower panel: Masks we adopt 
in Polarization. Blue region is for case "a". Blue plus light blue region is for case "b". Blue plus light 
blue and plus green region is for case "c" . Blue plus light blue plus green and plus orange region is for 
case "d". See also Table HB 

V. VARIANCE ANALYSIS 



In order to measure such a deviation from WMAP 9 fiducial model we do not take into 
account the estimator given in Eq. ([8]) since it is not sensitive to a lack of power because of its 
quadratic expression in terms of distance from the theoretical model Cj . Therefore we decide 
to naturally consider the expressions for the Variance of the CMB fluctuations in pixel space, 
which can be written as follows, see e.g. Appendix A of [31] 

(11) 



"max /np t -I \ 



i>2 



V 47r 



(Q^ 



E 

£>2 



2£+l 
Svr 



(cr 



+ C. 



bb\ 



(12) 



(f/2) = (Q2) 



(13) 



where (5T, is the fluctuation of the temperature map and Q and U are the Stokes parameters 
maps, see again e.g. [31] and references therein. 





FIG. 4. WMAP 9 year CMB spectra at low resolution. Red symbols for case "a", Black symbols for 
case "b" , Blue symbols for case "c" and Violet symbols for case "d" . See also Table |llj From top to 
bottom, TT, EE, TE, BB, TB and EB spectrum. For each spectrum we show £{£ + 1)Ci/2tt (given in 
/iK^) versus the multipole order (,. See also the text. 
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FIG. 5. Zoom of the TT spectrum given in Fig.|4]at the lowest muhipoles. See also Table |TT] 



A. Temperature analysis 



In Fig. p^ we show the variance {ST"^) computed through Eq. (11) with £max = 32 for all 
the cases of Table |TI} Each panel of Fig. [6] shows the histogram obtained with MC simulations 
and the vertical bar corresponds to the WMAP data. Fig. [6] shows that increasing the Galactic 
mask, the WMAP TT variance is more and more anomalous. The probability to obtain a 
smaller value than WMAP from the random simulations is 9.87%, 2.2%, 0.4% and 0.01% for 
the case "a", "b", "c" and "d" respectively 

In Fig. [7] we plot such a percentage for all the cases of Table [II| as a function of imax- This 
Figure makes clear that the behavior given in Fig. [6] for Imax = 32 is in fact a general one and 
it holds for every choice of imax from 2 to 32. In particular, for the largest mask we consider, 
i.e. case "d", the temperature variance is not consistent with ACDM model at 99.94% — 99.99% 
C.L. in the £max range [32,40]. 

The low amplitude of the TT variance that we measure for the WMAP 9 year data, is 
dominated by the lowest multipoles. In order to show this, we plot in Fig. [8] the probability for 
the TT variance to obtain a smaller value than the WMAP one for the cases "a" (solid line), 
"b" (long dashed line) , "c" (dashed line) and "d" (dotted line) of Table In] as a function of imin 
with imax = 32. Fig. [8] shows that excluding the quadrupole and the octupole in case "d" the 
anomaly disappears. 



Note that the more interesting cases "a" and "d" are analyzed with 10* random simulations whereas "b" and 
"c" with 10'^ random simulations. 
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FIG. 6. TT Variance computed with Eq. (11) with e,nax = 32 for the cases "a", "b", "c" and "d" from 



top to bottom and left to right respectively. See also Table H Each panel shows counts versus (ST^). 
The vertical bars correspond to WMAP 9 year data. Upper left and lower right panels are built with 
IC* random extractions whereas upper right and lower left with 10"^ random extractions. 
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FIG. 7. The probability for TT variance to obtain a smaller value than the WMAP one for the cases 
"a" (solid line) , "b" (long dashed line) , "c" (dashed line) and "d" (dotted line) of Table [ll] as a function 



of ' 



B. Impact of the smoothing procedure in the temperature analysis 



We check here the stability of our findings against the degradation procedure of the Temper- 
ature map, needed to perform the analysis at low resolution. Specifically we have degraded the 
WMAP 9 year ILC map to Ngi^^ = 16 performing no smoothing. Thanks to 10^ MC simulations 
we have re-computed Eq. (11) for cases "a" and "d" with (max = 32, see Fig. [dI The probability 



11 



80 - 



60 - 



40 



a 20 




/. 



U^'" 






6 8 10 12 

minimum multipole 



14 



16 



FIG. 8. The probability for TT variance to obtain a smaller value than the WMAP one for the cases 
"a" (solid line), "b" (long dashed line), "c" (dashed line) and "d" (dotted line) of Table [TT| as a function 

of train with ir,iax = 32. 



600 


kq85 


500 




^ " ^ FWHM=0 : 


400 




ri- |h : 


300 




1 


200 




1 - 


100 



-, ^- 


[T^r^„„_^ 



kq85+12deg 
FWHM=0 




1500 20(H) 



2500 3000 3500 4000 4500 5000 1500 2000 2500 



3500 4000 4500 5000 



FIG. 9. TT Variance computed with Eq. (11 1 with imax — 32 for the cases "a" and "d" from left to right 
respectively. See also Table pi] Each panel shows counts versus {^T"^). The vertical bars correspond to 
WMAP 9 year data of the ILC map pixellized at Nside — 16 without any smoothing. These panels are 
built with 10 random extractions. 



to obtain a smaller value than WMAP from the random simulations is 5.02% and 0.01% for the 
case "a" and "d" respectively. Therefore our results are proven to be stable. 

C. Polarization analysis 



We have repeated for polarization the same analysis we have performed for temperature, see 



Eqs. ( 12 ) , ( 13 ) , finding no evident anomaly. For sake of brevity we report in Fig. 10 the analogous 



plot for polarization, already given in Fig. [7] for temperature. Even if in a regime that is noise 
dominated already at £ > 10, it is interesting to notice that increasing the polarization mask, 
the percentages move towards less anomalous values in a monotonic way with respect to the 
sky fraction (at least for imax ~ 25). 
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FIG. 10. The probability for (Q^) to obtain a smaller value than the WMAP one for the cases "a" (solid 
line), "b" (long dashed line), "c" (dashed line) and "d" (dotted line) of Table III] as a function oi imax- 

VI. DISCUSSIONS AND CONCLUSIONS 



Using a QML estimator, we have provided new APS estimates of WMAP 9 year data at 
large angular scales and studied their dependence on Galactic masking. We analyzed all six 
CMB spectra. 

By using realistic MC simulations we have found that the lowest multipoles of the TT 
spectrum decrease in amplitude as the Galactic mask is increased. This indicates that the 
temperature anisotropics around the Galactic mask behave statistically differently from the 
anisotropics at high Galactic latitude. In particular, for the largest mask we consider, which 
covers ~ 54% of the sky, the temperature variance is not consistent with ACDM model at 
99.94% - 99.99% C.L. in the Imax range [32,40]. 

We have also checked that the degradation procedure to obtain the low resolution WMAP 9 
year ILC map is not impacting our findings. Specifically, we have degraded the WMAP 9 year 
ILC map to low resolution (corresponding to Nside = 16) performing no smoothing^ and we 
have re-evaluated for this map, Eq. ( |11[ ) in cases "a" and "d" . Our results are also stable against 
the smoothing procedure, and therefore there is no evident spurious effect that is introduced 
in this data manipulation step. Furthermore we have verified that the lowest multipoles are 
responsible of such a variance anomaly. 

Though based on different estimators, our findings are in agreement with the analysis per- 
formed in j3], for WMAP 5 and 7 data, and with [18 1, where WMAP 3 year data is considered. 



We remind that our results presented in Subsection VA are based on the WMAP 9 year ILC map smoothed 



with FWHM=9.1285 degrees, see Section III 
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We do not find evidence of a similar effect in the otlier polarized spectra. This is not surprising 
since they are noise dominated already at £ ~ 10. 

Recently Planck analysis of the variance for the temperature CMB anisotropics has ap- 
peared |26j . We aim at analyzing these data as well along the same lines as given here for the 
temperature case since the pattern of TT APS at large angular scales is close to the WMAP 9 
year data when the same mask is adopted ^22j. Note that the lack of power at large angular 
scales claimed in p2] is not the same as the one discussed here, since it involves different range 
of multipoles (see Fig. 39 of [22j). 

In conclusion, we speculate that this "low temperature variance anomaly" might be an 
indication of a non-standard phase of the cosmological evolution of the Inflaton before the slow 
roll phase as studied for example in [2], [6] and [5] for different pre-inflationary scenario^^' 
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